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We study numerically the effects of loop quantum gravity motivated corrections on massless 
scalar field collapse in Painleve-GuUstrand coordinates. Near criticality, the system exhibits Chop- 
tuik scaling with the added features of a mass gap and a new scaling relationship dependant upon 
the quantum length scale. The quantum corrected collapse exhibits a radiation-like phase which 
resolves the singularity: the black hole consists of a compact region of space-time bounded by a 
single, smooth trapping horizon. The "evaporation" is not complete but leaves behind an outward 
moving remnant. 
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Introduction. A successful quantum theory of gravity is expected to resolve the singularities 
present in general relativity. The two leading approaches to quantum gravity are string the- 
ory and loop quantum gravity (LQG). A method developed by Thiemann pp to regularize 
the gravitational potential in LQG has been successful in demonstrating singularity resolu- 
tion in mini-superspace models of black holes [2l, 13] and quantum cosmologies [3]. Moreover, 
it has long been conjectured that singularity resolution in black hole spacetimes is intimately 
related to Hawking radiation. Hayward [5] in particular has presented an elegant model of 
a complete non-singular black hole spacetime in which Hawking radiation at the apparent 
horizon prevents the formation of a static event horizon. The black hole is replaced by 
a compact region of trapped surfaces bounded by a single smooth trapping horizon. To 
the best of our knowledge, neither singularity resolution in LQG nor Hayward's radiative 
scenario has been explicitly realized in the context of a dynamical field theory. 

The purpose of the present letter is to provide such a dynamical realization by examining 
numerically the effect of a LQG motivated quantum correction to the gravitational potential 
in the context of spherically symmetric scalar field collapse. We work in Painleve-Gullstrand 
coordinates because they are regular across apparent horizons and allow for the evolution 
of trapped surfaces. The proposed quantum corrections are motivated by the Thiemann 
trick [Ij: instead of diverging at the origin, the gravitational potential is modified to vanish 
smoothly within some quantum length scale. The resulting numerical calculations produce 
a black hole spacetime that is non-singular and has a compact region of trapped surfaces 
instead of the usual event horizon. In addition, the quantum corrections lead to interesting 
new properties of Choptuik scaling [6], including a mass gap which was also observed 
recently by Husain [7|. 



Equations of Motion. We begin with a brief statement of the classical equations developed 
in [8]. The dynamical system that we study is that of a collapsing spherically symmetric, 
massless scalar field in four spacetime dimensions. 
The PG metric is given by 



ds^ = -aHe + [dr + ^^^acit j + rHn\ (1) 
where r'^dVt^ is the metric on a two sphere of radius r. The lapse function a and the Misner- 
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Sharp mass M. are determined by: 



Al' = l(.V)^ + ^V\/— -A'n*. (2) 



I2GM 



where prime denotes differentiation w.r.t. r. The dynamical equations for the scalar field tfj 
and its canonical conjugate 11^ are: 



, 2GM , 



o I rV + \\ — ^n. 



(5) 



These equations need to be supplemented by boundary conditions for M. and o . Without 
loss of generality we scale the time coordinate by choosing oit, 0) = 1. We also fix M.{t, 0) = 
which guarantees a flat metric in the neighbourhood of the origin. A non-zero value of M. 
at the origin signals the formation of a singularity. 

We now motivate quantum corrections that may arise due to the underlying spatial 
discreteness implied by LQG. We look to modify terms deriving from the potential which 
may diverge at the origin. These are easily identified as the 2GAA/r terms in Eqs.(|2]) to ([s]). 
The n^/r^ terms are kinematical and do not diverge. 

Polymer quantization [9J as used in LQG mini-superspace models replaces the momentum 
p conjugate to x by the finite translation operator = e*^^. The quantum version f/^ = e*^^ 
has an action given by: 

{U,^){x)=^{x + fi). (6) 

For the moment we consider x to span the entire real line. The Thiemann trick [Ij replaces 
the classical divergent function funreg = sgn{x)/^/\x\ by a commutator: 



Notice that freg is well defined at the origin but has a cusp at a; = /i. As in j7j, we instead 
use a smoothed function that has the same qualitative features as freg- it is nearly equal to 
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funreg for X S> /i, it has a positive sfope for x < /i and it is equal to zero at x = 0. 



f smooth (■^) 



sgn(xj 



\x\ 



1-e 



(8) 



Since the equations of motion depend on r > (rather than x) we revert to the radial 
coordinate. See Fig. [T] for a comparison of funreg, freg and f smooth in this domain. The 
quantum correction is simply the replacement 



2GM 



V2GMf smoother). 



(9) 



in each of (jSj-jSj). These corrections effectively add a repulsive component to the gravitational 
potential very near the origin, which is consistent with the expectation that the underlying 
discreteness implied by LQG at the microscopic level will give rise to a short range repulsion 
in the semiclassical limit. 




FIG. 1: (Color online.) A comparison of funreg{r), freg{r) and fsmoothi^) with = 0.1. 



Numerical Methods. Our code is essentially identical to that used in the recent investigation 
of Choptuik scaling in PG coordinates where further details can be found. One must 
first specify on a discrete spatial grid the initial data for the scalar field. Two forms were 
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considered: 

2 



tjj = Ar^ exp 



r — ro 
B 



(10) 



^ = ^tanh(^^-^j, (11) 

with n^(r, 0) = in all cases. A, B and tq are parameters which may be varied to study 
critical collapse. 

Equations (|2]) and (|3]) determine the mass and lapse functions which are inserted into the 
right hand side of the time evolution equations (|4]) and ([s]) in order to calculate the scalar 
field and its conjugate at the next time slice. These are then used as initial data for the 
subsequent evolution. This process is repeated until either a singularity forms or the matter 
disperses to leave behind a flat spacetime. 

In order to adequately resolve the small r structure we refined the r-spacing Ar(r) near 
the origin. The code destabilized if information moved over too many r-points in a single 
time step, so we used an adaptive time step At(t) refinement 

At(t) = MIN, I ^ Ar(r) | , (12) 

where ^ is the inverse of the local speed of an ingoing null geodesic and MIN,. denotes the 
minimum across the time slice. 



Quantum Corrected Choptuik Scaling. Initial data can be separated into two regions sepa- 
rated by a critical value, p^,, of some parameter p. If the ADM mass increases with increasing 
p, a black hole forms for p > p* with the inevitable singularity classically. If p < p*, no hori- 
zon forms and all matter disperses to infinity. If p ~ p*, the solution temporarily echoes the 
critical solution that describes a discretely self-similar zero mass black hole. The existence 
of this critical solution as an intermediate attractor gives rise to the universal mass scaling 
properties first observed by Choptuik [B]. 

Our quantum mass scaling results are given in Fig. [2j These curves in PG coordinates 
are non-sinusoidal with large amplitude cusps 0. The period and slope are in excellent 
agreement with previous studies of Choptuik scaling despite the unusual functional form. 
The new feature in the present case is a change in the critical parameter (which we refer to 
as Pg^) that is dependant upon the choice of quantum length scale /x. For sufficiently small 
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fi <^ Ar(0), the system is numerically identical to the classical case (see the = 10~^° line in 
Fig. |2]). For larger yU (though still small enough to be near criticality) the solutions exhibit 
critical scaling for some range p > Pg^ and come to an abrupt stop at p = pg^, below which 
no black holes form. This is the mass gap predicted by quantum gravity [2], and recently 
found by Husain [7] using a numerical model in null coordinates with similar corrections 
to those used here. It is interesting that in PG coordinates larger pg^ does not necessarily 
correspond to larger black hole mass due to the non-monotonic, large amplitude oscillations 
in the scaling relation. 



* n^lO'^-lO ■ |x=0.0005 n=O.OQl o n=0.005 □ n=0.01 * n=0.05 




I -12^ 

ln|A-A*| 

FIG. 2: (Color online.) Plots of the mass scaling behaviour for a variation in the parameter A in 
gaussian initial data on a mesh with Ar = 10"'* near the origin. 

Remarkably, we find near criticality a quantum power scaling relation: 

Pq* = kp{j) +P*^ (13) 

where kp is a family dependant constant, / is a length parameter required for proper units 
and (3 = 2.29 ± 0.04 is a universal constant whose value was verified by varying A and B/l 
in both forms of initial data. See Fig. [3] for an example. 
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FIG. 3: Variation in the parameter A in tanh initial data. 

Singularity Avoidance. Figure |4] provides spacetime diagrams for gaussian initial data far 
from criticality {p ^ Pq*) for the classical and quantum corrected potentials. Animations of 
the time evolution of the mass density corresponding to Figs. 4(a) and |4(b)' are available for 



viewing at http://theoryx5.uwinnipeg.ca/users/jziprick/ . In both cases inner and 



outer apparent horizons (where outward moving radial null rays are stationary) appear at 
some initial PG time and then separate as the inner horizon moves inward. In between the 
two horizons exist trapped surfaces on which all null rays move inward. For generic initial 
data the outer horizon moves outward until the last of the matter has fallen through it. 
Interestingly, a second pair of horizons forms for a period in both cases before annihilating 
to leave a single pair. 

In the classical case, the singularity forms when the inner horizon reaches the origin, at 
which point the code terminates. The outer horizon remains stationary once it has reached 
its maximum radius; it cannot decrease because the scalar field obeys the standard energy 
conditions. 



Figure 4(b) clearly exhibits singularity avoidance in the quantum case. As expected, 
the outer horizon rot reaches a maximum radius after all mass in the vicinity has passed 
through. As the inner horizon nears the quantum length scale it begins to slow down, 
eventually bouncing and moving outward without ever reaching the origin. Near the bounce, 
the mass density takes on negative values. This leads to mass loss and causes the outer 
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horizon to shrink until it meets the inner horizon in annihilation, leaving behind an outward 
moving remnant with strictly positive mass density. The amount of mass lost during the 
non-conservative "evaporation" phase seems numerically to be arbitrarily small for horizons 
forming in the quantum region which evaporate very soon after formation (for ^/r oh -max ~ 
1). For larger black holes {^i/roh. ,max ~ -^/^) observed mass losses of up to 80%, and 
higher values are likely possible. One can speculate that the quantum corrected collapse of 
macroscopic black holes would also leave microscopic remnants with radii of order fi. 

According to Einstein field equations, the energy density measured by an observer at 
constant r is 

M' 

P=^, (14) 



The negative mass density near the quantum bounce thus implies a violation of the null and 
dominant energy conditions inside the black hole. Such violations are also expected [11] in 
the presence of Hawking radiation. 
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FIG. 4: (Color online.) Spacetime diagrams. The (black) solid lines are null geodesies and (red) 
dashed lines indicate the trapping surface boundary. 

For a time after the outer horizon reaches its maximum and before it starts moving 
inward, the spacetime appears from the exterior to be that of a static black hole. The 
amount of time that the outer horizon remains stationary depends upon the distance the 
mass has to travel from r = Voh down to the quantum region near r ^. The shape of the 
mass distribution is also a factor since even for Toh ^ /x, the outer horizon may continue to 
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grow until the onset of evaporation if the infalhng mass is spread out over a distance of the 
same magnitude or greater than Voh- By adjusting the initial mass distribution and fi, the 
lifetime of the seemingly classical black hole can be made arbitrarily large or small. 

In the study by Hayward 0, a small length parameter was used to effectively preclude 
singularity formation. A Vaidya-like region with positive ingoing energy flux was matched to 
a stationary black hole spacetime to model black hole formation. After an arbitrary "time" 
given by the null coordinate v, a Vaidya-like region with negative ingoing energy flux was 
matched to the stationary black hole to model evaporation. The negative energy flux was 
balanced by outgoing Hawking radiation in order to conserve energy. This scenario bears 
a striking similarity to the black hole evolution produced by our numerical simulations, 
suggesting that the disappearing mass in our system is in some sense being radiated away, 
despite the fact that we have not included Hawking radiation directly. The key point in our 
analysis is that the quantum correction contains a repulsive core that prevents the inner 
horizon from reaching the origin, so that the singularity cannot form. The subtle (non-local) 
interplay between the microscopic behaviour of the potential and the dynamics of the outer 
horizon, accompanied by violations of the energy conditions, results in acausal behavior 
of the semiclassical system. Moreover, since there is a mass remnant left behind after 
"evaporation" , our semiclassical scenario has more than one feature relevant to the ultimate 
resolution of the information loss paradox. Although many questions remain, our anal- 
ysis provides tantalizing support for singularity avoidance in dynamical black hole formation. 
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